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The diagram technique for the one-band Hubbard model is formulated for the 
case of moderate to strong Hubbard repulsion. The expansion in powers of the 
hopping constant is expressed in terms of site cumulants of electron creation and 
annihilation operators. For Green's function an equation of the Larkin type is 
derived and solved in a one-loop approximation for the case of two dimensions 
and nearest-neighbor hopping. With decreasing the electron concentration in 
addition to the four bands observed at half-filling a narrow band arises near the 
Fermi level. The dispersion of the new band, its bandwidth and the variation 
with filling are close to those of the spin-polaron band in the t-J model. 
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1. Introduction 

Systems with strong electron correlations which, in particular, cuprate per- 
ovskites belong to are characterized by a Coulomb repulsion that is comparable 
to or larger than hopping constants. One of the simplest models for the description 
of such systems is the Hubbard model [ EJ El E] with the Hamiltonian 

H ^ ] ^nmfljjfffmff ~F ~ ^ ] ^-ncr^n,— cr; (1) 
nmcr ^ ncr 

where t nm is the hopping constants, the operator a no . creates an electron on the site 
n with the spin projection a = ±1, U is the on-site Coulomb repulsion and the 
electron number operator n na = a no .a n <7- in the case of strong electron correlations, 
U > |£ n m|, h is reasonable to use a perturbation expansion in powers of the hopping 
constants for investigating Hamiltonian ([T}. Apparently the first such expansion 
was considered in reference [Hj. The further development of this approach was given 
in references [ El El El E] where the diagram technique for Hubbard operators 
was developed and used for investigating the Mott transition, the magnetic phase 
diagram and the superconducting transition in the Hubbard model. 
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The rules of the diagram technique for Hubbard operators are rather intricate. 
Besides, these rules and the graphical representation of the expansion vary depend- 
ing on the choice of the operator precedence. The diagram technique suggested in 
references [ El El El E] is free from these defects. In this approach the power 
expansion for Green's function of electron operators a ncr and a]^ is considered and 
the terms of the expansion are expressed as site cumulants of these operators. Based 
on this diagram technique the equations of the Larkin type [E] for Green's function 
were derived [ El El E| • 

However, the application of this approach runs into problems. In particular, at 
half-filling the spectral weight obtained after a resummation of diagrams appears to 
be negative near frequencies ujd = ±U/2 [ 13J. This drawback is connected with di- 
vergencies in cumulants at these frequencies [E]- As can be seen from formulas given 
below, all higher-order cumulants have such divergencies at with sign-changing 
residues, which are expected to compensate the negative spectral weight in the entire 
series. On the other hand, at frequencies neighboring to u>d cumulants are regular. 
If a selected subset of diagrams is expected to give a correct estimate of the entire 
series for these frequencies the values at ujd can be corrected by an interpolation us- 
ing results for the regular regions. This procedure was applied in reference [ E] for 
the case of two dimensions, half-filling, nearest-neighbor hopping and with the use 
of the one-loop approximation. The spectrum was shown to consist of four bands. 
These band structure and the calculated shapes of the electron spectral function are 
close to those obtained in the Monte-Carlo [E1EIE], cluster perturbation [E] 
and the two-particle self-consistent [120] calculations. The four-band structure of the 
spectrum of the Hubbard model was also considered in reference [ |2T] ■ 

In the approach of reference [E] the mentioned four-band structure of the spec- 
trum at half-filling has its origin in the regions of large damping which separate the 
low- and high-frequency bands. In the major part of the Brillouin zone both these 
types of the bands owe their origin to the same terms of the irreducible part, i.e. 
to the same interaction processes. This brings up the questions: How this spectrum 
is changed with filling and how the quasiparticle peak, which determines the pho- 
toemission leading edge, arises in the spectrum? As expected, this peak differs in 
nature from other spectral features. In the present work it is shown that on certain 
deviation from half-filling an additional narrow band arises near the Fermi level on 
the background of the four above-mentioned bands. The dispersion of the new band, 
its bandwidth and the variation with filling are close to those of the spin-polaron 
band in the t-J model [ 1221123] . 

In the following section the perturbation expansion for the electron Green's func- 
tion is formulated in the form convenient for calculations and Larkin's equation is 
derived. In Sec. Ill the equations of the previous section are used for calculating the 
spectral function and the obtained results are discussed. Concluding remarks are 
presented in Sec. IV. 
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2. Diagram technique 

We consider Green's function 

G(nY, nr) = (TcVcr(T> n<T (r)>, (2) 

where the angular brackets denote the statistical averaging with the Hamiltonian 
TC — H — nJ2na n na, A* is the chemical potential, T is the time-ordering operator 
which arranges other operators from right to left in ascending order of times r, 
gw(t) = exp(7ir)a n(T exp(— TCr) and a n(J (r) = exp(7ir)a ncr exp(— TCt). Choosing 

H 0=-^Y1 n nan n - a - /i n ™ and H X = E ^nmaL a m CT (3) 
ncr ncr nmcr 

as the unperturbed Hamiltonian and the perturbation, respectively, and using the 
known expansion [|2I] for the evolution operator we get 

G(nV,nr) = E J ■ ■ ■ / ^...dr, E " E • • • ^ 

fc=0 ' ' ' nm^CTi n fe nj. cr fc 

x (Ta n / CT (r / )a nCT (r)a n / (7l (ri)a ni(Tl (ri) . . . an' fc cr fe (rA;)a nfc(Tfe (r fe ))oc, (4) 

where (3 = T~ l is the inverse temperature, the subscript "0" near the angular bracket 
indicates that the averaging and time dependencies of operators are determined with 
the Hamiltonian H . The subscript "c" indicates that terms which split into two or 
more disconnected averages have to be dropped out. 

The Hamiltonian Hq is diagonal in the site representation. Therefore the average 
in the right-hand side of equation (|3J) splits into averages belonging to separate lattice 
sites. To be nonzero these latter averages have to contain equal number of creation 
and annihilation operators. Let us consider the average which appears in the first 
order: (Ta n / CT (r')a nCT (r)a n ' . 1 (ri)a nicri (ri)) (for short here and below subscripts "0" 
and "c" are omitted). For this average to be nonvanishing operators have to belong 
either to the same site or to two different sites, 

(TUn'a (r) a na (r) a n / CTl (t[ ) a„ lCTl (n ) ) 
= (Ta nr7 (r) a na (r ) a nr7l (t[ ) a nr7l (r x ) ) 8 nn > 5^ 5 nni 

nicri nicri 

-(Ta ni0 -(r / )a ni0 - 1 (ri))(Ta no - 1 (r{)a ncr (r))5 n / ni 5 mi ' i (l - 5 nni ). 

The multiplier 1 — S nni ensures that n ^ ni in the two last terms (the case n = ni is 
taken into account by the first term in the right-hand side). After the rearrangement 
of the terms we find 

(Ta n / CT (r')ancr(i")a n ' 1 cr 1 (rOanicr 1 (ri)) = K 2 (T'a,Ta,T[a 1 ,Tia 1 )5 im >5 im ' i 5 IlIll 

+K 1 (t'ct, rcr)ifi(r{cri, nffi)^^^' - Ki(r'a, na^K^r^, ra)5 n / ni (5 nn ' , 
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where 

K\{j'o',to) = (Ta a (r')a a (T))8 a(7l , 

K 2 {t'o,to,t[o 1 ,t 1 o 1 ) = (Ta .(r / )o .(r)a fll (7^)a 01 (ri)) (5) 
-K^t'o, to~)Ki(t[o~i, noi) + K^t'o, no^K^Ox, to) 

are cumulants [ |2H] of the first and second order, respectively. All operators in 
cumulants (0) belong to the same lattice site. Due to the translational symmetry of 
Hq in equation Q the cumulants do not depend on the site index which is therefore 
omitted in equation (jSJ). 

Averages which appear in 
higher-order terms of equati- 
on (JIJ can be transformed in the 
same manner. The average in 
the k-th order term which con- 
tains k + 1 creation and annihi- 
lation operators is represented 
by the sum of all possible prod- 
ucts of cumulants with the sum 
of orders equal to k + 1. All 
possible distributions of opera- 
tors between the cumulants in 
these products have to be taken 
into account. The sign of a term 
in this sum is determined by 
the number of permutations of 
fermion operators, which bring 
the sequence of operators in the 
initial average to that in the 
term. Actually the above state- 
ments determine rules of the diagram technique. Additionally we have to take into 
account the presence of topologically equivalent terms - terms which differ only by 
permutation of operators -£/i(rj) in equation (jlj). Since these terms are equal, in the 
expansion only one of them can be taken into account with the prefactor v = j/k\ 
where j is the number of topologically equivalent terms. Following reference [ [H] in 
diagrams we denote a cumulant by a circle and the hopping constant t nn / by a line 
directed from n' to n. The external operators a^^ij') and a nCT (r) are denoted by 
directed lines leaving from and entering into a cumulant. The order of a cumulant 
is equal to a number of incoming or outgoing lines. Summations and integrations 
over the internal indices rij, n^, Oi and are carried out. Since site indices of oper- 
ators included in a cumulant coincide, some site summations disappear. Also some 
summations over cij get lost, because in any cumulant spin indices of creation and 
annihilation operators have to match. Taking into account the multiplier (— l) fc in 
equation (@J), the sign of the diagram is equal to (—1)' where / is the number of loops 
formed by hopping lines. Figure ^ demonstrates connected diagrams of the first four 
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Figure 1. Diagrams of the first four orders of 
expansion Q. 
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orders of the power expansion (j3J) with their signs and prefactors. Here the thick 
line with arrow in the left-hand side of the equation is the total Green's function. 
Notice that if we set t nn = 0, contributions of the diagrams (b), (d)-(f), (i)-(n), 
(p), (t), and (u) vanish. However, below a renormalized hopping parameter will be 
introduced which is nonzero for coinciding site indices and therefore the mentioned 
diagrams are retained in figure E 

All diagrams can be separated into two categories - those which can be divided 
into two parts by cutting some hopping line and those which cannot be divided in 
this way [ |Hl EI] • These latter diagrams are referred to as irreducible diagrams. In 
figure[T]the diagrams (c), (e), (f), (h), (j), (k), (n), (p), and (r)-(v) belong to the first 
category, while the others are from the second category. The former diagrams can be 
constructed from the irreducible diagrams by connecting them with hopping lines. 
Notice that a prefactor v of some composite diagram is the product of prefactors 
of its irreducible parts [cf., e.g., the irreducible diagrams (a), (d) and the composed 
diagrams (j ) , (k) in figure ^ . If we denote the sum of all irreducible diagrams as 
fT(nV, nr), the equation for Green's function can be written as 

G(nV, nr) = K(n'r', nr) + W" dr^nV, n 1 r 1 )t nin; G(n / 1 r 1 , nr). (6) 

/ ^ 

This equation has the form of the Larkin equation [I15j. 

The partial summation can be carried out in the hopping lines of cumulants by 
inserting the irreducible diagrams into these lines. In doing so the hopping constant 
t nn ' in the respective formulas is substituted by 

0(nr, nV) = t nn ,5(r - r') + £ t n< f dr 1 K(n / 1 r, n 1 r 1 )0(n 1 r 1 , nV). (7) 

/ J 

For the diagram (b) in figure E inserting irreducible diagrams into the hopping line 
leads to the diagrams (g), (m), and (q). 

Due to the translation and time invariance of the problem the quantities in 
equations © and (J7|) depend only on differences of site indices and times. The use 
of the Fourier transformation 

G(k, iut) = T e ik(n '- n) f dr'e iWi(T '- r) G(nV, nr), 

Jo 

where lo\ = (21 + l)7rT is the Matsubara frequency, simplifies significantly these 
equations: 



l-t k A"(k,iw,) 



e(k> z ) = * k , =t k + tlG(kM)- (9) 

1 - t k K(k, iui) 

In the approximation used below the total collection of irreducible diagrams 
K(k, iui) is substituted by the sum of the two simplest diagrams (a) and (b) ap- 
pearing in the first two orders of the perturbation theory. Due to the form of the 
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latter diagram this approximation is referred to as the one-loop approximation. In 
the diagram (b) the hopping line is renormalized in accordance with equation 0. 
Thus, 

K(ko{) = Kx(iui) - K 2 (iLJia,iui 1 a 1 ,iui 1 a 1 )^J2tlGik,iLJi 1 ), (10) 

Zicti k 

where Ki(iui) and 

= JjjJ drdrd^drie ; i 1 1 K 2 [r o.to.t^^txOx) 

= p5i + i ljVJrl > x K 2 (iui<j, ic^ai, iio h ai) 

are the Fourier transforms of cumulants (jSJ), N is the number of sites and we set 
Sk^k — 0. Notice that in this approximation K does not depend on momentum. 

Now we need to calculate the cumulants in equation (JTUJ). To do this it is conve- 
nient to introduce the Hubbard operators X% = \in)(jn\ where |zn) are eigenvectors 
of site Hamiltonians forming H , equation (JHJ). For each site there are four such states: 
the empty state |0n) with the energy E Q — 0, the two degenerate singly occupied 
states | an) with the energy E\ = —fi and the doubly occupied state |2n) with the 
energy E 2 = U — 2/i. The Hubbard operators are connected by the relations 

a na = X^ + aX n °> 2 , aL = X? + aX 2 n ^ (11) 

with the creation and annihilation operators. The commutation relations for the 
Hubbard operators are easily derived from their definition. Using equation (jllj) the 
first cumulant in equation (jSJ) can be computed straightforwardly: 

1 /e-^+e-^o e-^+e~^\ 

Kiim) = — — + — , 12 

Z V iw, - E a0 iuji - E 2a J 

where Z = e~P E ° + 2e~^ Ea + e'^ 2 is the partition function, E tj = Ei - Ej. As 
indicated in [0 EH EB], if K{k,iui) is approximated by this cumulant the result 
corresponds to the Hubbard-I approximation 

To calculate K 2 it is convenient to use Wick's theorem for Hubbard operators [ 

(TX ai (n) . . .X ai (r,)X Q (r)X Qi+1 (r m ) . . . X an (r n )) 

= j2(-l) Pk 9 a (r - T k )(TX ai { Tl ) . . . [X ak ,X a ] ± (r k ) . . .X an (r n )), (13) 
fc=i 

where the averaging and time dependencies of operators are determined by the 
Hamiltonian H , a is the index combining the state and site indices of the Hubbard 
operator. If X a is a fermion operator (X 0cr , X u2 and their conjugates), P k is the 
number of permutation with other fermion operators which is necessary to transfer 
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the operator X a from its position in the left-hand side of equation (fTT?|) to the position 
in the right-hand side. In this case 



e&yT ( 

e m 3 ZjZ 1 1 e 11 



;1, r>0, 
;< J, r < 0, 



(14) 



where z and j are the state indices of X a . If X a is a boson operator (X 00 , X 22 , X au ' , 
X 02 , and X 20 ), P fc = and 



1, r>0, 
e^y, r < 0. 



(15) 



In equation ([T3]h [X afc ,X w ]± denotes an anticommutator when both operators are 
of fermion type and a commutator in other cases. 



The substitution of equation (JTT|) in (Ta CT (r')a a (r) 



(ti)) in X 2 , equa- 



tion (|5|). leads to six nonvanishing averages of Hubbard operators (such averages 
are nonzero if the numbers of the operators X° CT and X CT ° coincide in them, and 
the same is true for the pair X°" 2 and X 2cr ). Applying Wick's theorem (fT??J) to these 
averages the number of operators in them is sequentially decreased until only time- 
independent operators are left. For H in equation (J3J) these are X 00 , X ffCT ', and X 22 . 
Their averages are easily calculated. As a result, after some algebra we find 

^2K 2 (iuJia,iuJi 1 a 1 ,iuJi 1 ai) = 

-Z~ l U {e- (3Eo g 0a (iuj l )g 0r7 (iuji 1 )g 2(iuJi + iu h ) [g a(^i) + 9oA iuJ h 



+e 



+e 



-PE 2 



g<r2(m)g*2(iuj h )go2(iui + iu h 



g^ui) + g< 7 2{m 1/ 



g 0a [m)ga2 [g a fait ) 
+9oa(m 1 )g r72 (iu; h )(g^(iuJi) + gl 2 (m))] } 
-Z^ 2 U 2 p5 lh (e-^ Eo+E2) + 2e-^ Eo+E ^ + 3e" 2/3Bl + 2e~^ +E ^ 

X 90a( iuJ l)9l2( iuJ l) 

+Z 2 U 2 l3{2e~ l3{Eo+E2 ' ) + e -P( E o+ E ±) _|_ e -/3(Ei+E 2 ) 

x (iw/ ) g a 2 (m ) goa fah) g<,2 (iwh ) , 



(16) 



where gijiiuji) = (iui + -Eij) -1 is the Fourier transform of functions ([Hj) and (|15|). 

Equations (fT2Jl and (fTB^) can be significantly simplified for the case of principal 
interest U T. In this case if satisfies the conditions 



A < ^ < U — A, 



(17) 



A ^> T, the exponent exp(— /3-Ei) is much larger than exp(— /3£?o) an d exp(— /?i? 2 ). 
By passing to real frequencies we can ascertain that terms in (jlfjjl with the two latter 
multipliers contain the same peculiarities as other terms. Therefore terms with these 
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multipliers can be omitted in the equations and we get 

K ^ luJ i) = T- ; — vr- ; m> 

(iui + fj,){icui + fi-U) 

Y J K 2 {iL0i(7,\u h (T 1 ,\u h (j l ) = --Ug 0a (iuJi)g a 2(^i) 9oA iuJ h) +^2(^/1) 

1 r 12 3 



;is) 



-^Ug Qa {\u h )g a2 {\u h ) 



g Qu [iui) - g^KW 



Let us turn to real frequencies by substituting iui with z = u + irj where 7] is a 
small positive constant which affords an artificial broadening. Results given in the 
next section were calculated with G(k, z) in the right-hand side of equation (fTU|) 
taken from the Hubbard-I approximation. As mentioned, this Green's function is 
obtained if K(ku) in equation (jSJ) is approximated by K\ from equation (|18|) which 



gives 




U 2 + tl Z ~ £ 2,k 



(19) 



ei,k = \ (u + t k + y/lP + tl) - /x, e 2M = l(u + t±- yJU 2 + tl) - fi. 

Below the two-dimensional square lattice is considered. It is supposed that only 
the hopping constants between nearest neighbor sites t are nonzero which gives 
tk = 2t[cos(k x ) + cos(fc^)] where the intersite distance is taken as the unit of length. 
Due to the electron-hole symmetry of Hamiltonian (0) the consideration can be 
restricted to the range of the chemical potentials \i < ¥. 



3. Spectral function 

Figure El demonstrates QK(u) calculated with the use of equations (fTO|). (|TH|) 
and (|19|) . The change to real frequencies carried out in the previous section converts 
the Matsubara function (J2J) into the retarded Green's function [121] • It is an analytic 
function in the upper half-plane which requires that QK(uj) be negative. As seen 
from figure |2J this condition is violated at uJd = — A* an d U — fx. This difficulty 
of the considered approximation was indicated in reference [ Hi] . The problem is 
connected with divergencies at these frequencies introduced by functions go a (u)) and 
g a 2{ UJ ) m the above formulas. As can be seen from the procedure of calculating 
the cumulants in the previous section, these functions and divergencies with sign- 
changing residues appear in all orders of the perturbation expansion (@J). In the entire 
series the divergencies are expected to compensate each other so that the resulting 
QK(uj) is negative everywhere. However, in the considered subset of terms such 
compensation does not occur. Nevertheless, as seen from figure [21 at frequencies 
neighboring to uj^ cumulants are regular and, if the used subset of diagrams is 
expected to give a correct estimate of the entire series for these frequencies, the value 
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i 1 rn 1 r 




M oVU 



Figure 2. The imaginary part of K(uj) calculated using equations (|1U|) . (|T%|) 

and JEI) for a 100x100 lattice, i = -U/8 and T = O.OOlf/ (the dashed lines), 
(a) /i = 0.5C7, 7? = O.Olt/. (b) fi = 0.1U, rj = 0.02C7. The solid lines show the 
corrected QK(uj) (see text). 



of ^KIuj) at the singular frequencies can be reconstructed using an interpolation 
and its values in the regular region. Examples of such interpolation are given in 
figure El 

As seen from figure^, at half-filling, = ~, Q\fT(u;) has two broad minima. With 
decreasing the chemical potential from this value these minima shift with respect 
to the Fermi level without a noticeable change of their shapes until the Fermi level 
enters one of the minima at fi ~ 0.17U. As this takes place, two new sharp minima 
arise near frequencies —fi and U — fi on the background of the above-mentioned 
broad minima. The appearance of the broad features in figure 121 is connected with 
the third term in K 2 in equation ()18|) . while the sharp minima are related to the 
second term in this formula. Its contribution to K(lv), equation (J1U|) . grows rapidly 
when the Fermi level enters the broad minimum. 

The function K(z) has to be analytic in the upper half-plane also and therefore 
its real part can be calculated from its imaginary part using the Kramers-Kronig 
relations. We use this way with the interpolated QK(u) to avoid the influence of the 
divergencies on $tK(uj). However, the use of the interpolation overrates somewhat 
values of 1 3-^(^)1 which leads to the overestimation of the tails in the real part. To 
correct this defect the interpolated K(u) is scaled so that in the far tails its real 
part coincides with the values obtained from equation (fTU|) . 

The spectral function 

AQsu) = --QGQsu) = ^ r ^f 9 ( ^ r — - . 19 (20) 

obtained in this way for momenta along the symmetry lines of the square Brillouin 
zone is shown in figure El The shapes of the spectral function in figure are nearly 
the same as at half-filling - with decreasing [i from 0.5U to approximately 0.17U 
these curves shift with respect to the Fermi level without perceptible changes in 
their shapes. As indicated in reference [E], four bands can be distinguished in these 
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12- 



(0,0) 



(7C.0) 



(71, 7l) 



(0,0) 



(d) 











(0,0) 



(jt,0) 



{n,n) 



(0,0) 



Figure 3. Panels (a) and (c): the spectral function A(k.u) calculated for momenta 
along the symmetry lines of the square Brillouin zone in a 100x100 lattice for 
t = -[7/8, T = 0.00117, r] = 0.02J7, /i = 0.2J7 (a) and n = 0.05C7 (c). The 
dispersions of maxima in the panels (a) and (c) are shown in the panels (b) and 
(d), respectively. Here darker areas correspond to larger intensities of maxima. 
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spectra. For parameters of figure these bands are located near frequencies — 4|t|, 
\t\, 4|t|, and 9|i| (see figure Efc>). For the major part of the Brillouin zone the peaks 
forming the bands arise at frequencies which satisfy the equation 1 — t]JR:K(uj) = 
and fall into the region of a small damping l^i^u;)! [see equation (|20|) ]. As seen from 
figure 121 such regions of small damping are located between and on the outside of 
the two broad minima in QK(u). This is the reason of the existence of the four well 
separated bands - two of them are located between the minima of QK(u), while two 
others are on the outside of these minima. Broader maxima of A(ku) for momenta 
near the boundary of the magnetic Brillouin zone are of different nature - since is 
small for such momenta, the resonant denominator in equation EH does not vanish 
and the shape of the spectral functions is determined by QK(u) in the numerator 
of this formula. 

More substantial changes in v4(kco>) occur for fi < 0.17U when the Fermi level 
enters one of the broad minima in QK(u). As seen from figure Efc, in addition to the 
mentioned four bands there appear sharp dispersive features near u = — \x and U — ji 
for momenta in the vicinity of the boundary of the magnetic Brillouin zone. It is clear 
that these changes in the spectral function are connected with the changes in ^sK{uS) 
shown in figure El The peaks near — /i are more intensive than those near U — fx and 
are located in the nearest vicinity of the Fermi level. For \i w 0.17C/ the width of 
the band formed by the former peaks is comparable to the superexchange constant 
J = At 2 /U of the effective Heisenberg model which describes magnetic excitations in 
the limit U ^> \t\. This indicates obviously the participation of the spin excitations in 
the formation of these band states. The bandwidth decreases with further reduction 
of the electron concentration. As this takes place, the peak intensities first grow and 
then saturate. The maximum energies of the band are located near the boundary 
of the magnetic Brillouin zone. For parameters of figure Efc for these momenta the 
band touches the Fermi level and the corresponding peaks disappear above it (see 
figure Eli) . 

These properties of the band resembles those of the spin-polaron band in the t-J 
model. This latter band is also located near the Fermi level, has the similar dispersion 
and the bandwidth, which decreases with the reduction of the electron concentration 
(with the increase of the hole doping counted from half-filling) [EH|- However, there 
is one essential difference in the behavior of these bands in the Hubbard and t-J 
models. In the former model the narrow band near the Fermi level appears at a 
certain deficiency of electrons, while the above-mentioned four bands exist in the 
entire considered range of electron concentrations. In the t-J model the situation 
is opposite - the spin-polaron band exist in the wide range of hole concentrations 
< x < 0.17, while the wider band - an analog of the four-band structure - starts 
to form at x « 0.06 [123]. 

The value of the chemical potential for which the Fermi level enters the mini- 
mum of QK(uj) and the narrow band begins to form depends on the ratio t/U. For 
example, for t = —U/A this happens at fi w 0.27U. 

There is some difficulty in the comparison of our results with the data of Monte- 
Carlo calculations and cluster theories. It is connected with the fact mentioned in 
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reference [Q3h the one-loop approximation overestimates the spectral weights of the 
two internal bands of the four-band structure near the momenta (0,0) and (tt,tv). 
As a consequence the electron concentration calculated from the formula 



appears to be considerably underestimated. For example, for the parameters of fig- 
ures OK andEfc such estimated (n) equals to 0.87 and 0.68, respectively. From the 
comparison of the dispersions in figures and Eli with the results of Monte-Carlo 
calculations (figure 9 in [HE]) it can be concluded that these concentrations have 
to be approximately 0.95 and 0.85, respectively. With this in mind we find that 
the spectral functions and dispersions in figure El are close to the Monte-Carlo spec- 
tra (cf. with figures 10 and 11 in [EH])- In these latter spectra for some deviation 
from half-filling a weakly dispersive feature is also observed near the Fermi level. 
However, this feature has low intensity and is lost at the foot of a more inten- 
sive maximum on approaching the boundary of the magnetic Brillouin zone. These 
differences may be connected with the comparative high temperature T = 0.33 |t| 
used in the Monte-Carlo simulations. A similar weakly dispersive feature near the 
Fermi level was obtained also in quantum cluster theories, however for much smaller 
deviations from half- filling (cf. with figure 2c in [|2Zj and figure 32 in [I28j). 

4. Conclusion 

The considered diagram technique is very promising for a generalization to many- 
band Hubbard models for which energy parameters of the one-site parts of the Hamil- 
tonians exceed or at least are comparable to the intersite parameters. The expansion 
in powers of these latter parameters can be expressed in terms of cumulants in the 
same manner as discussed above. Now there are distinct cumulants which belong to 
different site states characterized by dissimilar parameters of the repulsion and the 
level energy. These cumulants are described by formulas similar to equations (|T2~j) 
and (jlbp. For example, for the Emery model [ 12*9*] which describes oxygen 2p a and 
copper 3d x 2_ y 2 orbitals of Cu-0 planes in high-T c superconductors there are two 
types of cumulants corresponding to these states. Diagrams of the lowest orders, 
e.g., for Green's function on copper sites resemble those shown in figure ^ where 
"oxygen" cumulants are included in hopping lines. Equations of the type of (JBJ) and 
partial summations similar to equation Q can be derived also in this case. 

In summary, the diagram technique for the one-band Hubbard model was for- 
mulated for the case of moderate to strong Hubbard repulsion. The expansion in 
powers of the hopping constant is expressed in terms of site cumulants of electron 
creation and annihilation operators. For Green's function the equation of the Larkin 
type was derived and solved for the case of two dimensions and nearest-neighbor 
hopping. With decreasing the electron concentration in addition to the four bands 
observed at half-filling a narrow band arises near the Fermi level for momenta near 
the boundary of the magnetic Brillouin zone. On the occurrence the width of the 
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band is comparable to the super exchange constant J = At 2 /U which indicates the 
participation of the spin excitations in the band formation. The bandwidth decreases 
with decreasing the electron concentration. The maximum energies of the band are 
located near the boundary of the magnetic Brillouin zone. For some deviation from 
half-filling in these points the band touch the Fermi level. By these properties the 
band resembles the spin-polaron band of the t-J model. 



References 

1. Gutzwiller M.C., Phys. Rev. Lett., 1963, 10, 159. 

2. Hubbard J., Proc. R. Soc. London, Ser. A, 1963, 276, 238. 

3. Kanamori J., Prog. Theor. Phys., 1963, 30, 275. 

4. Hubbard J., Proc. R. Soc. London, Ser. A, 1966, 296, 82. 

5. Westwanski B., Pawlikovski A., Phys. Lett. A, 1973, 43, 201. 

6. Slobodyan P.M., Stasyuk I. V., Teor. Mat. Fiz., 1974, 19, 423 (in Russian) [Theor. 
Mat. Phys., 1974, 19, 616]. 

7. Zaitsev R.O., ZhETF, 1976, 70, 1100 (in Russian) [Sov. Phys. JETP, 1976, 43, 574]. 

8. Izyumov Yu.A., Skryabin Yu.N., Statistical Mechanics of Magnetically Ordered Sys- 
tems. Consultants Bureau, New York, 1988. 

9. Ovchinnikov S.G., Valkov V.V., Hubbard operators in the theory of strongly correlated 
electrons, Imperial College Press, London, 2004. 

10. Vladimir M.I., Moskalenko V.A., Teor. Mat. Fiz., 1990, 82, 428 (in Russian) [Theor. 
Mat. Phys., 1990, 82, 301]; Vakaru S.I., Vladimir M.I., Moskalenko V.A., Teor. Mat. 
Fiz., 1990, 85, 248 (in Russian) [Theor. Mat. Phys., 1990, 85, 1185]. 

11. Metzner W., Phys. Rev. B, 1991, 43, 8549. 

12. Moskalenko V.A., Entel P., Digor D.F., Phys. Rev. B, 1999, 59, 619. 

13. Pairault S., Senechal P., Tremblay A.-M.S., Eur . Phys. J. B, 2000, 16, 85. 

14. Sherman A., Preprint arXiv:cond-mat/0601331 , 2006. 

15. Larkin A.I., ZhETF, 1959, 37, 264 (in Russian) [Sov. Phys. JETP, 1960, 37, 186]. 

16. Moreo A., Haas S., Sandvik A.W., Dagotto E., Phys. Rev. B, 1995, 51, 12045. 

17. Preuss R., Hanke W., von der Linden W., Phys. Rev. Lett., 1995, 75, 1344. 

18. Grober C, Eder R., Hanke W., Phys. Rev. B, 2000, 62, 4336. 

19. Dahnken C. et al., Phys. Rev. B, 2004, 70, 245110. 

20. Tremblay A.-M.S., Kyung B., Senechal D., Preprint [ar"Xiv:cond-mat/0511334[ 2005. 

21. Shvaika A.M., Phys. Rev. B, 2000, 62, 2358; Stasyuk I.V., Shvaika A.M., Ukrainian 
J. Phys., 2002, 47, 975. 

22. Plakida N.M, Oudovenko V.S., Yushankhai V.Yu, Phys. Rev. B, 1994, 50, 6431. 

23. Sherman A., Schreiber M., Phys. Rev. B, 1994, 50, 12887; Pseudogaps and mag- 
netic properties of the two-dimensional t-J model. - In: Studies of High Tempera- 
ture Superconductors, vol. 27, p. 163. Nova Science Publishers, New York, 1999, ed. 
Narlikar A.V.; Eur. Phys. J. B, 2003, 32, 203. 

24. Abrikosov A. A., Gor'kov L.P., Dzyaloshinskii I.E., Methods of Quantum Field Theory 
in Statistical Physics. Pergamon Press, New York, 1965. 

25. Kubo R., J. Phys. Soc. Jpn., 1962, 17, 1100. 

26. Moukouri S. et al., Phys. Rev. B, 2000, 61, 7887. 

27. Kyung B. et al., Preprint |arXiv:cond-mat/0502565 , 2005. 



13 



A. Sherman 



28. Maier T. et al., Rev. Modern Phys., 2005, 77, 1027. 

29. Emery V.J., Phys. Rev. Lett., 1987, 58, 2794. 



14 



